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Use of the Sunyaev-Zeldovich effect as a precise cosmological probe neces- 
sitates a realistic assessment of all possible contributions to Comptonization of 
the cosmic microwave background in clusters of galaxies. We have calculated 
I the additional intensity change due to various possible populations of energetic 

electrons that have been proposed in order to account for measurements of intra- 
cluster radio, nonthermal X-ray and (possibly also) EUV emission. Our properly 
normalized estimates of (the highly model dependent value of) the predicted in- 
tensity change due to these electrons is well below ~ 6% and ~ 35% of the usual 
Sunyaev-Zeldovich effect due to electrons in the hot gas in Coma and A2199, 
respectively. These levels constitute high upper limits since they are based on 
energetic electron populations whose energy densities are comparable to those 
of the thermal gas. The main impact of nonthermal Comptonization is a shift 
of the crossover frequency (where the thermal effect vanishes) to higher values. 
Such a shift would have important consequences for our ability to measure clus- 
ter peculiar velocities from the kinematic component of the Sunyaev-Zeldovich 
effect. 
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1. Introduction 

The Sunyaev-Zeldovich (S-Z) effect is a small intensity change that results from Comp- 
ton scattering of the cosmic microwave background (CMB) radiation by electrons in the hot 
gas in clusters of galaxies (Zeldovich & Sunyaev 1969, Sunyaev & Zeldovich 1972). The 
effect constitutes a unique cluster and cosmological probe (for reviews, see Rephaeli 1995a, 
Birkinshaw 1999), whose great potential is beginning to be realized in recent years follow- 
ing the major observational progress in obtaining sensitive images of the effect by (mostly) 
interferometric arrays (Jones et al. 1993, Carlstrom et al. 2000, Udomprasert, Mason, & 
Readhead 2000, and the review by Carlstrom et al. 2001) and the increasingly more precise 
values of the Hubble constant {H ) that have been deduced from S-Z and X-ray measure- 
ments. For example, a fit to data from 33 cluster distances yields Hq = 58 km s^ 1 Mpc -1 (in 
a flat cosmological model), with direct observational errors of ±5% (Carlstrom et al. 2001), 
but with a much larger level of systematic uncertainties. Among others, the latter include 
simplified modeling of the properties of the hot intracluster (IC) gas, and cluster morphology. 
These are the main limitations to the use of the S-Z effect as a precise cosmological probe, 
and as such give further motivation for in-depth studies of the cluster environment. 

Energetic non-thermal (NT) electrons whose pressure is not negligible compared to 
the thermal gas pressure constitute an aspect of IC phenomena with possibly appreciable 
ramifications for precision S-Z measurements. The presence of significant energetic electron 
populations in many clusters has been known from measurements of diffuse IC radio emission, 
and recently also from NT X-ray emission in a few clusters (Rephaeli, Gruber & Blanco 1999, 
hereafter RGB, Fusco-Femiano et al. 1999, Kaastra et al. 1999, Gruber & Rephaeli 2002). 
NT electrons produce an additional degree of Comptonization which amounts to a small 
intensity change (AI nt ) that must be accounted for, particularly in the measurement of H Q 
from the thermal component, and peculiar cluster velocities from the kinematic component 
of the S-Z effect. Relativistic generalizations (Rephaeli 1995b, Challinor & Lasenby 1998, 
Sazonov & Sunyaev 1998) of the original non-relativistic calculations of Sunyaev & Zeldovich 
(1972, 1980) have now been performed to a sufficiently high level of accuracy, including also 
terms of order r 2 (Nozawa, Itoh, & Kohyama 1998, Itoh et al. 2000, Shimon & Rephaeli 
2002), where r is the Thomson optical depth of the cluster. In calculating the effect of 
multiple scatterings the finite size of the cluster has to be explicitly accounted for; this has 
been done in the Monte-Carlo simulations of Molnar & Birkinshaw (1999). The relativistic 
treatment provides the theoretical basis for calculation also of AI nt . 

CMB Comptonization by NT electrons was first assessed for conditions in lobes (McK- 
innon, Owen, & Eilek 1991) and cocoons (Yamada & Sugiyama 1999) of radio galaxies, with 
its possible use to measure their electron pressure. The effect of such electrons in clusters 
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was recently considered in some more detail (Blasi & Colafrancesco 1999, Ensslin & Kaiser 
2000, Blasi, Olinto, & Stebbins 2000). Clearly, the higher the electron pressure, the higher is 
the degree of Comptonization induced by the electrons, and because models of NT electrons 
vary greatly in energy density, estimates of their impact on the CMB range from a negligible 
level, to a very substantial fraction of the thermal S-Z effect due to the hot gas. It is quite 
essential to study NT electron populations in order to realistically determine the spectral 
and spatial profiles of AI nt in clusters in which extended radio emission has been measured. 
This may lead to ways by which their impact on S-Z work can be minimized, and eventually 
the Comptonization induced by NT electrons could perhaps even be used as a diagnostic 
tool of these electrons. 

In this paper we perform an exact calculation of the degree of Comptonization predicted 
in a range of NT electron models that have been proposed to explain radio, EUV, and hard 
X-ray emission in clusters. Our estimates of AI nt in Coma and A2199 are based on recent 
results on the electron populations in these clusters, and the main features of the resulting 
spectral change are contrasted with those of the thermal component of the S-Z effect. 

2. Non-thermal Electron Populations 

The main evidence for relativistic electrons and magnetic fields in clusters comes from 
measurements of extended IC regions of radio emission (in the frequency range ~ 0.04 — 4 
GHz) with spectral indices and luminosities in the range ~ 1 — 2, and 10 40 ' 5 - 10 42 erg/s 
(H = 50 km s _1 Mpc -1 ), respectively. Many of the radio emitting regions detected in over 
30 clusters (Giovannini, Tordi, & Feretti 1999, Giovannini & Feretti 2000) are central, with 
sizes in the range ~ 1 — 3 Mpc. The energy range of the emitting electrons depends on 
the value of the mean, volume-averaged magnetic field, a quantity which is not known very 
well, but is likely to be in the range ~ 0.1 — 1 /iG (for a recent review, see Rephaeli 2001). 
The range of electron energies implied from these measurements is roughly 1 — 100 GeV, 
but electrons with energies both below and above this range are also expected on theoretical 
grounds. 

Relativistic electrons produce X-ray emission in a wide spectral range by Compton scat- 
tering off the CMB (Rephaeli 1979) and by NT bremsstrahlung. This emission has quite 
possibly been measured already in Coma (Rephaeli, Gruber & Blanco 1999, Fusco-Femiano 
et al. 1999), A2199 (Kaastra et al. 1999), A2256 (Molendi, De Grandi, & Fusco-Femiano 
2000), and A2319 (Gruber & Rephaeli 2002) by the RXTE and BeppoSAX satellites. These 
were spectral measurements; the PCA, HEXTE (both aboard RXTE) and PDS (aboard Bep- 
poSAX) experiments do not have the adequate spatial resolution. Thus, we have no spatial 
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information on NT electrons from X-ray measurements, and only rudimentary knowledge of 
the morphology of the radio emission (which, however, involves also the unknown spatial 
distribution of the magnetic field). Therefore, in order to avoid the need for introducing un- 
known parameters, we characterize electron populations in terms of the spectral distribution 
of their total number, ignoring their spatial profiles in the central ~ 1 Mpc region where 
they mostly reside. 

Of the various proposed NT IC electron models (see, e.g. , Sarazin 1999) we focus here on 
those that have been contrasted with actual observational data, with the electron distribution 
appropriately normalized by the determination of both their total number and the power-law 
index, q. Some of the proposed electron models were motivated by the presumption that 
NT X-ray emission could also originate from NT bremsstrahlung by energetic - though not 
necessarily highly relativistic - electrons, and by an attempt to account also for observed 
EUV emission from a few clusters. This emission is said to be NT (e.g. , Sarazin & Lieu 
1998, Bowyer et al. 1999), possibly by a population of low energy electrons. The full electron 
distribution could be described either as a sum of two separate parts (thermal plus NT), or 
by a 'super Maxwellian', consisting of a truncated Maxwellian with a power law tail (Blasi, 
Olinto, & Stebbins 2000) 'sawn' together at a given energy (e.g. , ~ 3kT e , where T e is the 
electron [and gas] temperature). We consider here four specific models whose parameters in 
Coma and A2199 have been determined from radio, EUV, and NT X-ray measurements. 

The simplest model for the electron momentum distribution is a power-law over a suffi- 
ciently wide range so as to explain both the observed radio and possibly Compton-produced 
X-ray emission. We express the electron spectrum in terms of the (differential) number, 
N(p), per unit dimensionless momentum, p = Pj, where 7 and f3 are the Lorentz factor and 
dimensionless velocity, (3 = v/c, respectively, 

N(p) = Ap~ q , pi < p < p 2 

_ N l2 (q - 1) 

- „:<»->-„-<»->' (1) 

The limiting momenta p\ and p 2 are the lower and higher values of p that correspond to 
the observed radio frequency range, and N 12 is the total number of electrons with energies 
in this specific interval, [pi,p2]- Details of such a model were worked out long ago (Rephaeli 
1977, 1979) and will not be repeated here. Suffice it to say that the full momentum range 
is substantially uncertain, especially its low end which is of particular relevance to our 
discussion here. 



A simple way to obtain the lower energy extension of the above distribution is to assume 
that steady state is attained, whereby electrons are continually accelerated to compensate 
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for radiative energy losses - defined here in terms of dp/ dt - and a population of lower energy 
electrons is built up. The latter can be determined by taking into account the dominant rate 
of energy loss at energies well below ~ 150 MeV, electronic (Coulomb) excitations, given by 
(Rephaeli 1979) 



b PP = 1.2 x 10 _1 V 



1.0 + 



In (Vl+p 2 /n e ) 



75 



(2) 



where n e is the (thermal) electron number density. At higher energies Compton-synchrotron 
losses dominate; these occur at a combined rate (e.g. , Rephaeli 1979) 



b cs = 1.37 x 10- 20 [(1 + zf + 9.5 x 10- 2 (5/buG) 2 ] s _1 , 



(3) 



where z is the cluster redshift, and B is the mean, volume-averaged value of the magnetic 
field. The steady state solution is then (Rephaeli 1979) 

Ab P 



Nj(p) = 



Hp) 

Ab cs 
b{p) 



(1+P 2 ) 
(1+P?)" 



(g-2) 
2 



(9-2) 
2 



(pi < P < P2) 
(p < Pi) , 



(4) 



where b = dp/dt = b ee + b cs (1 +p 2 )- Clearly, the observed radio frequency range implies 
that pi 1. (Note that another possible energy loss mechanism is scattering by Alfven 
waves, which the electrons themselves can excite if their spatial distribution is somewhat 
anisotropic. Because of the substantial uncertainty in estimating this loss mechanism [e.g. , 
Rephaeli 1979] it will not be considered here.) Parameters of this model were determined 
from RXTE measurements of the Coma cluster, assuming the emission is from Compton 
scattering of relativistic electrons off the CMB (Rephaeli, Gruber, & Blanco 1999), and 
from BeppoSAX measurements of A2199 (Kaastra et al. 1999). This appropriately extended 
electron population which has a power-law form at low and high energies, but with indices 
whose values differ by unity, is our first model. The parameters of all the models considered 
here are given in Table 1. 

Two other forms for N(p) are discussed by Sarazin & Kempner (2000). In the first model 
electrons are accelerated to relativistic energies by shell-type supernova remnants (Baring 
et al. 1999). When the back reaction of the accelerating electrons on the shock structure 
is taken into account, then this presumably leads to a modified power-law form at small 
momenta, resulting in the following distribution, 



N n (p) = 



Npp- 

i+p! 



1+ ^ 



(5) 
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Cluster 


Model 


Model details 


N 




N nt 




N nt/ N ih 

lit 1 ill 




Coma 




















I 


Extended power law 


2.275 x 


10 72 


1.09 x 


10 65 


4.70 x 10" 


-7 




II 


Shock accelerated 


4.270 x 


10 67 


2.98 x 


10 68 


1.28 x 10" 


-3 




III 


Cooling electrons 


5.238 x 


10 62 


3.13 x 


10 65 


1.35 x 10" 


-6 




IV 


Simple power law q=3.68 


8.970 x 


10 68 


9.27 x 


1Q 69 


0.04 




A2199 




















I 


Extended power law 


2.306 x 


1Q 68 


1.66 x 


1Q 64 


2.10 x 10" 


-7 




II 


Shock accelerated 


4.240 x 


10 67 


6.22 x 


10 68 


7.90 x 10" 


-3 




III 


Cooling electrons 


5.238 x 


10 62 


3.13 x 


10 65 


3.95 x 10" 


-6 




IV 


Power law q=3.33 


5.420 x 


10 68 


6.47 x 


1Q 69 


0.082 





Table 1: Models for NT electrons in the Coma and A2199 clusters. 

We assume p c = 1 in normalizing the model. Explicit account of Coulomb losses yields the 
third distribution 

AT ( \ 2N °P 2 (*\ 

N is derived from comparison of the models with observations and assuming that the hard 
X-ray emission mechanism is NT bremsstrahlung. The electron distribution in Equation (6) 
has an upper momentum cutoff at p — 300. 

Sarazin & Kempner (2000) also propose a simple power law model which extends to 
lower energies with no change in the index. This is the fourth (IV) model explored here. The 
observationally deduced parameters for the above four models are listed in Table 1; in these 
models, the lower momentum cutoff is pi = + 30) 2 — 1, where = kT e /mc 2 . 

3. Comptonization 

A detailed description of the Comptonized spectrum resulting from scattering of the 
CMB by a thermal, non-relativistic population of electrons was given by Sunyaev & Zeldovich 
(1972, 1980). A more accurate treatment of the process requires relativistic generalization 
due to the rapid motion of electrons in the hot (kT e < 15 keV) IC gas, as has been shown 
explicitly by Rephaeli (1995a). A fully relativistic treatment (but still in the Thomson limit, 
and with electron recoil and induced scattering safely ignored) is obviously required in order 
to calculate the Comptonized spectrum resulting from scattering of the radiation by a NT 
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energetic electron population. We use the approach adopted in the latter paper (see also 
Rephaeli & Yankovitch 1997). 

The probability of scattering of an incoming photon moving originally in the direction 
specified by /i = cos# > to the direction fi' = cos9' , is (Chandrasekhar 1950) 



/(// ,//o) = g 



(7) 



where the subscript denotes quantities in the electron rest frame. The logarithmic frequency 
shift in the scattering is 



s = In (v' I v) = In 



(8) 



where (3 is the dimensionless electron velocity in the CMB frame. The probability for scat- 
tering is conveniently written in terms of the variables f3 and t = e s (Wright 1979), 



V(s,P) = 



e S 



where /ii and fj, 2 are given by 



^PJ^ (l + M) 2 



e " 3(1 / ) - 1 S <0 

-1 s > 



(9) 



(10) 



1^2 



1 s < 

"-(y- 1 s >o 



Integration over /i in Equation (9) yields 

311 -t 



V{t,p) 



+ 



32p 6 t 

3(l + t) 
8p 5 



[1 + (10 + 8/ + V)t + t 2 ] 



3 + 3p 2 + p A 3 + 2p 2 



2p 



(2sinh" 1 (p) - | In 



where 



(11) 



(12) 



ln(i) | < 2sinh" 1 (p) 



(13) 



Compton scattering results in a frequency shift, x — > a;', where a; is the dimensionless 
frequency a; = hu/kT, and T is the CMB temperature. The corresponding change of the 
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photon occupation number through a pathlength along the cluster, An(x), is obtained by 
integrations over the scattering probability and electron momentum distributions, 

An <*> = T £ C dpdt ^ v (t - p) fe^r - ^t) ' (14) 

where r is the optical depth of the cluster to Compton scattering, and N tot is the total number 
of electrons in the cluster. The measured quantity is the change in intensity, AI(x)/i = 
x 3 An(x), where i = 2(kT) 3 /(hc) 2 . We use the latter equation to calculate the additional 
intensity change induced by energetic electrons, AI nt (x), in the electron models described 
in the previous section, with parameter values sampling the observationally deduced ranges, 
as specified in Table 1 for all the four electron models. 

If Comptonization by NT electrons is not taken into account it could affect the observa- 
tionally deduced value of the crossover frequency, xo, defined as the frequency at which the 
purely thermal effect vanishes. The exact value of xo has practical significance, since obser- 
vations near this frequency are optimal for the determination of cluster peculiar velocities 
(Rephaeli & Lahav 1991) from measurements of the kinematic component of the S-Z effect. 
We have calculated the modified value of the frequency where now the sum of the thermal 
and NT intensity changes vanishes. For this we used the analytic expression of Nozawa, 
Itoh, & Kohyama(1998) for the thermal and kinematic components of the S-Z effect, and 
computed xq using the measured values of the temperatures in Coma and A2199. 



4. Results 

The impact of each of the above four NT electron models depends very much on the 
total number of these electrons as determined by radio, NT EUV and X-ray emission. The 
thermal gas parameters are based on X-ray measurements. For Coma, recent XMM mea- 
surements yield kT e = 8.2 ± 0.4 keV for the gas temperature (Arnaud et al. 2001), and 
(since quantitative results for the gas density profile from either XMM or Chandra are not 
yet available) we have taken the ROSAT deduced values (Mohr, Mathiesen, & Evrard 1999), 
n e0 ~ (3.12. ± 0.04) x 10~ 3 cm" 3 , r c = 0.366 Mpc, and (5 n ~ 0.705, for the central electron 
density, core radius, and index in the expression for (the commonly used) density profile, 
n e (r) = n e o(l + r 2 /r 2 )~ 3/3n / 2 , respectively. The corresponding values in A2199 are kT e = 4.7 
keV, central electron density n e0 = 7.21 x 10~ 3 cm~ 3 , r c = 0.196 Mpc, and (3 n = 0.78 (Kaastra 
et al. 1999). These values are based on H = 50 km s _1 Mpc -1 . 

Clearly, most of the intensity change AI nt (x) is due to the more numerous low energy 
electrons. More specifically, a typical CMB photon has dimensionless frequency x = 2.701, 
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and the scattered photon frequency is on average 

- = 1 + \v 2 - (15) 
x 3 

The observationally relevant range of values of x is at most x < 20, and since the (the 
68% likelihood) interval of rms values of x\ Ax' ~ x^2p 2 /3 + 46p 4 /45, is very wide for 
p > 1, it follows that electrons with momenta much larger than at most a few tens are 
largely irrelevant for boosting photons to frequencies such that x < 20. (Obviously, all 
electrons contribute to the value of the Thomson optical depth and the decrement on the 
Rayleigh- Jeans side.) 

We have calculated AI nt (x)/io and total electron energies in the extended power-law 
model described by equation (4), using the observationally deduced model parameters in 
Coma and A2199. We emphasize that the basic model is a power law at high energies 
(> 1 GeV) with the index determined from radio measurements. The distribution is then 
appropriately extended to lower energies (by taking electronic excitations losses into account) 
resulting in a change in the value of the index at low energies. The power law indices 
and overall normalizations were taken from Rephaeli, Gruber & Blanco (1999) for Coma, 
and Kaastra et al. (1999) for A2199. Based on these parameters, the deduced values of 
AI n t(x)/io are very small in comparison with the magnitude of the thermal S-Z effect in 
these clusters as predicted from the measured values of the gas temperature and density. 
In Coma the S-Z effect was actually measured by Herbig et al. (1995) using the OVRO 
5.5 m telescope, and more recently by De Petris et al. (2002) using the MITO telescope. 
In the former paper the temperature change in the center of Coma was reported to be 
AT = —505 ± 92 fiK at 32 GHz (x ~ 0.56), a value which is consistent with that predicted 
based on the X-ray deduced parameters. A somewhat lower value was deduced by De Petris 
et al. (2002). Our calculated values of AI nt {x) /i$ due to NT electrons are smaller than 
0.001% of the predicted thermal effect in Coma and A2199. Clearly, the implied shift in the 
value of the crossover frequency is also negligible. The very small impact of NT electrons in 
this model is not surprising given the relatively low energy content of these electrons, ~ 0.8% 
of the thermal electron energy in Coma and A2199. 

Next we calculated the degree of Comptonization in the shock accelerated and cooling 
electrons models (II and III) described in equations (5) and (6). We have normalized these 
models to match the observed EUV emission, at a luminosity level of roughly 5 x 10 43 erg s _1 
over the band 0.065-0.245 keV in both Coma and A2199. In the shock accelerated electron 
model the calculated levels of AI nt (x) are < 0.5% and < 3% of the corresponding intensity 
change due to thermal electrons, and the shift in the value of the crossover frequency is 
also quite small. Calculated values of the energy in NT electrons in this model are ~3% 
and 13%-16% (for p c = 0.3 — 1.0) of the energy in thermal electrons in Coma and A2199, 
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respectively. In the cooling electron model, the degree of NT Comptonization is negligible, 
less than 0.003% of the magnitude of the thermal S-Z effect in both Coma and A2199. This 
is mainly due to the small number of NT electrons in this model, ~ 0.1% of that in the 
shock accelerated electron model. The NT electron energy constitutes 0.9% and 4.2% of the 
electron thermal energy in the two clusters. Our calculated quantities are listed in Table 2. 
The first column is the cluster name (or Abell number), the second is the model number, 
third is the ratio of (total) energy in NT to energy in thermal electrons. In the fourth column 
we list the net energy deposition rate (in keV/Gyr) by NT electrons, and in the fifth the 
value of the crossover frequency; values of the error in the peculiar velocity and the ratio 
AI nt (x) / AI are listed in last two columns. 



Cluster 


Model 


Ent/ E t h 


dE/dt 




Av > 


AI nt /AI 






% 


(keV/Gyr) 




(km/sec) 


% 


Coma 
















Thermal 






3.9000 








I 


0.79 


-0.06 


3.9000 




1.4 x 10" 4 




II 


3.18 


10.80 


3.9020 




0.41 




III 


0.90 


-0.04 


3.9000 




4.7 x 10~ 4 




IV 


21.59 


296.86 


3.9216 


-110 


6.8 


A2199 
















Thermal 






3.8700 








I 


0.74 


-0.15 


3.8701 




1.2 x 10" 4 




II 


15.87 


58.87 


3.8803 




2.97 




III 


4.23 


-0.06 


3.8702 




2.5 x 10 -3 




IV 


64.17 


626.39 


3.9385 


-210 


34.5 



Table 2: The impact of the four energetic electron models (see the text for definitions of the 
listed quantities). 

The much higher NT electron energy content in the single power law (model IV) sug- 
gested by Sarazin & Kempner (2000), ~ 20% - 22% (for q = 2.92 - 3.68) in Coma, and 
~ 64% — 184% (for q = 2.2 — 3.33) in A2199, result in substantial degrees of additional 
Comptonization. Using values of the parameters as deduced by Sarazin & Kempner (2000), 
listed in Table 1, we calculate the change of intensity, AI nt (x)/io, shown in Figure 1, and 
listed in Table 2. It is clear from this figure that AI nt (x) can reach an appreciable fraction 
of AI(x). For example, near the peak of the Comptonized Planckian, at x ~ 6.5, AI nt /AI 
is < 13% and < 35% for Coma and A2199, respectively. The implied shift in the value of 
the crossover frequency in this model, ~ 2 GHz in Coma, and > 4 GHz in A2199, would also 
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be observationally relevant. Such a systematic shift would introduce an error in the deduced 
value of the peculiar velocity. The implied error (which does not depend on the velocity, if 
the very small quadratic dependence of AI(x) on the velocity is ignored) amounts to ~ 110 
km/s for Coma (q = 2.92 - 3.68), and ~ 210 km/s (q = 3.33) and 320 km/s (q = 2.2) for 
A2199 (values lower than 100 km/s are not shown in Table 2). 




Fig. 1. — The spectrum of the additional Comptonization due to energetic electrons in the 
models proposed by Sarazin & Kempner (2000) for the Coma (left panel) and A2199 (right 
panel). 

Finally, we have calculated the heating rate of the hot gas by the NT electrons through 
Coulomb interactions (Rephaeli 1979, Rephaeli & Silk 1995) in the four models considered 
in this paper. Energy deposition rates into the gas were estimated by taking a mean gas 
density in the central 1 Mpc region of Coma and A2199. Subtraction of the observed cooling 
rate due to thermal bremsstrahlung (X-ray) emission then yields the net rate of temperature 
change. Values of the latter quantity (in keV/Gyr) are listed in Table 2. 



5. Discussion 



We have briefly discussed models for NT energetic electrons and their predicted effect 
on the CMB based on parameters that were deduced from radio, EUV, and X-ray mea- 
surements of Coma and A2199. Our immediate objective has been the exact calculation of 
the additional degree of Comptonization due to these electrons, and its possibly relevant 
observational consequences. Only electron models that have direct observational motivation 
were considered. Since all these models have basically a power-law form, it is clear that their 
impact on the CMB is largely due to the more abundant lower energy electrons (except for 
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the cooling electrons model). However, the interpretation of the EUV and X-ray emissions 
as NT bremsstrahlung and Compton scattering are not secure, so there is considerable un- 
certainty in quantifying electron densities. Moreover, since low energy electrons lose energy 
mostly by non-radiative coupling to the thermal gas, the uncertainty in the estimation of 
their density is particularly substantial. Only the presence of relativistic electrons is well 
established from measurements of extended regions of radio emission in many clusters (in- 
cluding Coma, but not A2199). Clearly, radio measurements yield the electron density only 
if we have an independent estimate of the mean, volume-averaged value of the magnetic field 
in the emitting region. Estimates of the mean IC field are few and uncertain (Newman, 
Newman & Rephaeli 2002). 

Due to the lack of detailed information on IC energetic electron populations, relevant 
theoretical considerations are of particular interest. First, since it is virtually always the case 
that particle acceleration mechanisms tap only a small fraction of the energy at the source, 
energy density in NT electrons is likely to be only a small fraction of the gas (which is a 
significant mass component of clusters) thermal energy density. Realistically, therefore, the 
thermal energy density in clusters is expected to set an absolute upper limit to the energy 
density in NT electrons, a limit which is very unlikely to be reached. Second, the energetic 
electron population is linked to the gas mainly by Coulomb coupling that results in energy 
transfer from energetic electrons (and protons) to the gas. This coupling and the resulting 
heating set stringent constraints on the attainment of steady state in general, and the density 
of low energy NT electron models in particular. 

We have calculated the ratio of total energy in NT electrons to that in thermal electrons 
in the four energetic electron models considered in this paper; for the first three models the 
calculated values (in Table 2) are low. For the power law model of Sarazin & Kempner 
(2000), the corresponding values are very high (~ 22% in Coma, and > 64% in A2199) and 
quite unrealistic. Such large energy contents would also imply a high rate of energy transfer 
to the gas, and heating at a rate higher than the rate of cooling by emission of thermal 
radiation. This without even considering the additional heating by IC energetic protons 
(whose energy density in the Galaxy is higher than that of relativistic electrons). 

In fact, the single power law model of Sarazin & Kempner (2000) is not only unappealing 
from an energetic point of view, but is also questionable on a more basic physical ground: 
As we have stated in Section 2 - based on the original work of Rephaeli (1979) - correct 
extrapolation of the relativistic electron energy spectrum to low (< 100 MeV) energies 
must account for electronic Coulomb excitations, the dominant energy loss mechanism at 
such energies for typical gas densities in the central regions of clusters. The very weak 
energy dependence of this mechanism (in contrast with the quadratic dependence of the 
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Compton-synchrotron energy loss rate) flattens the spectrum at low energies, resulting in a 
much smaller number of electrons than would have been predicted otherwise. We conclude 
that this model is non-viable; therefore, the high degree of implied Comptonization and 
appreciable shift in the value of Xq that are predicted in this model are at best very high 
upper limits to the impact of NT electrons on measurements of the S-Z effect. 

Our basic result in this paper is that exact calculation of the impact of realistically 
normalized models of NT electrons in Coma and A2199 yields levels of the degree of Comp- 
tonization by such electrons that are only a negligible fraction of the corresponding S-Z effect 
due to the hot IC gas. It is clear from our discussion that this result is generally valid since 
the energy density in NT electrons is not likely to be a significant fraction of the thermal 
energy density. In clusters with a significant NT electron population the added Comptoniza- 
tion due to these electrons would clearly constitute a source of confusion in the analysis of 
S-Z measurements. Since there are little or no differences in the spectral CMB signatures 
of thermal and NT electrons, high spatial resolution X-ray, S-Z and radio measurements of 
these clusters would be needed in order to minimize this confusion when (as expected) the 
spatial profiles of the two electron populations are found to be detectably different. 

Useful comments made by the referee on an earlier version of the paper are gratefully 
acknowledged. This research has been supported by the Israeli Science Foundation grant 
729\00-02 at Tel Aviv University. 
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